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摘要 : 其 长 基线 干涉 测量 ( VLBI ) 通过 相位 校准 ( PCAL ) 信号 来 校准 链 路 时 延 。 但 现 有 系统 的 PCAL 信号 只 能 
获取 链 路 时 延 的 相对 变化 。 而 面 对 UT1 精准 测量 、 台 站 钟 差 补偿 、 维 护 诊断 和 提升 未 来 深 空 探测 精度 等 需求 ， 天 
线 链 路 的 绝对 时 延 定 标 至 关 重要 。 在 天 线 链 路 中 ， 由 于 存在 变频 器 等 器 件 ， 绝 对 时 延 测量 较为 困难 。 本 文 利用 梳 状 
谱 相位 随 频率 线性 变化 的 特点 ,提出 了 一 种 不 需要 参考 变频 器 ,采用 PCAL 信号 定 标 绝对 链 路 群 时 延 的 方法 。 同 时 ， 
本 文 设计 了 一 款 用 于 台 站 快速 诊断 和 绝对 时 延 定 标的 轻 量 级 PCAL 提取 软件 ,该 软件 目前 已 在 实际 调试 工作 中 应 用 。 
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在 甚 长 基线 干涉 测量 ( VLBI ) 中 ,由 线 缆 、 变 频 器 等 导致 的 链 路 时 延 是 重要 误差 之 一 ,因此 需要 对 其 进行 校准 , 
才能 获得 真正 的 几何 时 延 。 现 有 校准 技术 为 在 前 端 接 收 机 处 注入 等 间隔 的 点 频 信 号 ， 即 相位 校准 ( PCAL ) 信号 ， 
相关 处 理 时 将 PCAL 信和 号 的 相位 提取 出 来 ， 通 过 计算 得 到 群 时 延 变化 曲线 ， 即 可 获取 准确 的 链 路 时 延 变化 ， 从 而 进 
行 校准 [9。 在 一 般 VLBI 观测 中 ，, 仅 关注 链 路 时 延 的 相对 变化 ， 绝 对 链 路 时 延 采用 差分 观测 消除 。 但 面向 未 来 ， 我 
国 深 空 探测 的 目标 是 要 向 更 深 、 更 远 的 天 体 进 行 研究 ， 需 要 更 精准 的 测 角 精度 ， 定 标 绝对 链 路 时 延 对 消除 误差 、 提 
升 精度 有 一 定 帮 助 。 在 UT1 精准 测量 、 台 站 钟 差 补偿 和 维护 诊断 中 ， 也 需要 对 链 路 的 绝对 时 延 进 行 定 标 018.191。 


对 于 抛物 面 天 线 的 结构 ， 按照 信号 传输 的 路 径 可 以 分 为 反射 面 光 程 段 、 馈 源 段 和 线 缆 段 。 其 中 光 程 段 时 延 可 通 
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过 几何 计算 得 到 ,其 因 俯仰 重力 和 温度 引起 的 变化 ， 可 通过 摄影 测量 2 和 微波 全 息 测量 进行 定 标 2224。 馈 源 段 主要 
由 喇叭 和 馈 源 网 络 组 成 ， 其 中 喇叭 段 时 延 较 稳 定 ， 可 忽略 不 计 2526]。 对 于 馈 源 网 络 和 线 缆 段 ， 比 较 复 杂 ， 因 此 需要 
PCAL 信号 进行 定 标 、 监 测 和 校准 。 

对 于 线 线 、 放 大 器 和 滤波 器 等 器 件 ， 可 以 用 矢量 网 络 分 析 仪 定 标 绝 时 延 。 对 于 变频 器 ,目前 的 方法 主要 为 引入 
参考 变频 器 , 或 者 通过 梳 状 谱 校准 标量 混 频 器 ,使 之 蔡 代 参考 变频 器 [12-171。 本 文 提出 一 种 用 PCAL 信号 定 标 绝对 链 
路 时 延 的 方法 ， 该 方法 利用 梳 状 谱 相 位 随 频 率 线性 变化 的 特点 ， 不 需要 参考 变频 器 ,同时 用 改变 频率 间隔 的 方法 解 
决 了 模糊 度 问 题 ， 可 直接 对 整个 天 线 链 路 绝对 时 延 进行 定 标 ， 用 于 满足 台 站 建设 和 维护 诊断 中 的 需求 。 

同时 ,目前 VLBI 台 站 所 用 PCAL 提取 软件 都 是 基于 大 型 相关 机 设计 ， 有 较 强 的 并 行 计算 能 力 和 批 处 理 能 力 ， 
在 处 理 多 台 站 PCAL 信号 时 展现 出 强大 的 性 能 [8-111。 但 面向 台 站 快速 诊断 和 本 文 提 到 的 绝对 时 延 定 标 等 方面 的 实际 
需求 ， 尤 其 是 提高 诊断 和 定 标 效率 的 需求 ， 轻 量化 软件 是 非常 必要 的 。 因 此 ， 本文 就 根据 该 需求 ， 设计 了 一 款 轻 量 


化 的 PCAL 提取 软件 ， 该 软件 目前 已 在 实际 调试 工作 中 应 用 。 


2 基本 原理 
PCAL 发 生 器 为 一 个 窄 脉冲 形成 器 ( 或 称 梳 状 频谱 发 生 器 ), 该 器 件 的 频率 基准 由 毛 原 子 钟 提供 ， 最 后 形成 宽 
度 为 毫 微 秒 量 级 的 脉冲 串 ,在 频 域 表现 为 等 频率 间隔 的 梳 状 点 频 信号 , 即 PCAL 信和 号。 根据 PCAL 发 生 器 的 原理 ， 
信号 的 时 域 表 达 式 为 : 
x(t) =X ô(t-nT) 
其 中 ,T 为 脉冲 周期 。 
通过 傅 里 叶 变 换 ， 可 得 : 


na SY (o 


根据 傅 里 叶 时 移 特 性 ， 经 时 延 rz， 时 域 和 频 域 分 别 变 为 : 
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x(D) = > 6(t- c nT) 
根据 上 式 ，PCAL 信号 中 ， 梳 状 谱 各 谱 线 的 相位 随 频率 呈 线 性 变化 。 因 此 ， 不 同 频 率 的 PCAL 位 点 可 以 得 到 相 
同 的 相 频 斜率 。 本 文 利 用 此 特性 ， 让 不 同 频率 的 PCAL 信号 互相 关 ， 求解 群 时 延 。 
如 图 1(a) 所 示 , PCAL 信号 发 生 器 在 天 线 的 馈 源 之 后 ， 前 级 放大 之 前 注入 ( PCAL1 ) , 与 观测 到 的 射电 信号 通 
过 同样 的 路 径 进 入 记录 器 ， 本文 将 该 链 路 称 为 待 测 链 路 ， 与 射电 信号 有 同样 的 链 路 时 延 。 同 时 ,在 记录 器 处 还 有 一 
个 PCAL 发 生 器 ( PCAL2 ), 该 PCAL 信号 经 防 混 翅 带 通 滤波 器 滤波 后 直接 进入 记录 器 ,本文 将 该 链 路 称 为 参考 链 
路 。 两 个 PCAL 发 生 器 采用 同一 本 振 参 考 ， 并 用 稳 相 传输 系统 锁 住 ， 保 证 初始 相位 一 致 。 同 时 ， 还 需 搭 建 一 校准 链 


”路 ， 如 图 1(b) 所 示 ， 用 于 消除 参考 链 路 的 固定 时 延 差 。 
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Fig. 1 Diagram of PCAL link 


根据 以 上 描述 ,由 于 PCAL 信号 的 频率 基准 为 所 原子 钟 ， 各 频率 分 量 均 为 谐 波 ， 相 位 是 相干 的 。 因 此 ， 在 注入 


处 ， 各 频 点 相位 为 : 
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(i) = ifop (0) 
其 中 ， 万 为 基 频 ，p(0) 为 基 频 初始 相位 
PCAL 信号 经 过 传输 线 缆 、 变 频 器 等 路 径 ,产生 r 的 时 延 ， 此 时 各 频 点 相位 为 : 
PÙ) = ifog(0) — 2nf, ` iT 
对 相 邻 PCAL 位 点 求 相 位 差 ， 可 得 : 
Aq = 9(i—1)— (i) = 2nfs1 + fop(0) 

分 别 对 待 测 链 路 与 参考 链 路 ， 校 准 链 路 与 参考 链 路 作 相关 处 理 ， 再 对 相 邻 PCAL 位 点 求 相 位 差 ， 可 得 : 

Apans 一 Aprer = 27fo(Tins 一 Ter) 十 无 [pi(0) -9z(0)] 

AQcai 一 Aprer = 27fo(Tear — Tref) + fole (0) — 692(0)] 


其 中 ,rins 为 待 测 链 路 的 时 延 , 即 真 正 的 绝对 链 路 时 延 。rrer 和 recu 分 别 为 参考 链 路 和 校准 链 路 的 时 延 ,其 中 rcuw 可 
相对 容易 地 用 矢量 网 络 分 析 仪 单独 定 标 。%i(0) — bz(0) 为 PCAL1 和 PCAL2 之 间 的 基 频 初始 相位 差 ， 由 于 采用 稳 
相传 输 系 统 锁定 ， 该 差 值 相对 稳定 。 


因此 将 两 式 相 减 、 移 项 便 可 得 到 绝对 链 路 时 延 : 


Agpins 一 Apecar 


Tins — 十 Tcal 
2n f, 
o 


在 实际 测量 中 ， 还 需 考虑 模糊 度 问 题 。 由 于 Aqp 的 范围 为 (n,n)， 因 此, RER EAA OC, E 


时 真实 的 时 延 可 以 表示 为 : 


"d 
T = tmesure T F 
fo 
其 中 : M 为 含 任意 自然 数 ， 为 PCAL 产生 的 模糊 度 
下 表 给 出 了 不 同 频率 间隔 下 ，PCAL 产生 的 模糊 度 : 
表 1 不 同 频率 间隔 的 模糊 度 


Table 1 Ambiguity of different frequency intervals 


Interval Ambiguity Measurement range Cable length 
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500KHz 2us (~1,1)us 480m 


1MHz 1us (—500,500)ns 240m 
2MHz 500ns (-250,250)ns 120m 
5MHz 200ns (7100,100)ns 48m 


其 中 : 群 速度 取 0.8c， 即 2.4 x 108m/s， 并 假设 纯 线 线 时 延 ， 无 变频 器 等 其 余 时 延 的 情形 
考虑 到 实际 台 站 的 线 缆 长 度 一 般 不 会 超过 100m , 总 链 路 时 延 一 般 不 会 超过 1us ,最 大 不 会 超过 2us , 因此 一 


般 情 况 下 1MHz 间隔 已 经 够 用 ， 可 保留 500KHz 设置 用 于 极端 情况 下 的 测量 。 


3 精度 评估 
如 图 2 所 示 , 横 轴 ( 实 轴 ) 与 纵 轴 ( 虚 轴 ) 分 别 表示 信号 向 量 在 实 部 与 虚 部 的 振幅 分 量 ， 信 号 向 量 与 实 轴 夹 角 
则 为 相位 。5 为 信号 向 量 ，N 为 噪声 向 量 ， 骨 为 合 向 量 ，8 向 量 与 邮 向 量 的 相位 差 即 为 误差 。 当 PCAL 信和 号 的 信 噪 比 


为 SNR， 积 分 时 间 为 7 时 ,单个 PCAL 位 点 相位 精度 为 : 


1 
n= 
?  J2TSNR 


根据 上 节 描 述 ， 时 延 本 质 为 相位 对 频率 的 斜率 。 实 际 测量 操作 为 ， 对 不 同 频率 处 带 有 误差 的 相位 作 线性 拟 合 。 


因此 ，, 拟 合 误差 即 为 时 延误 差 ， 有 : 


Qi = 2nfit +b 
Jacobian 系数 矩阵 为 : 
27 亡 1 
2nf. 1 
观测 误差 矩阵 为 : 
X= obl 


经 矩阵 变换 ， 可 得 时 延 精度 为 : 
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2 单个 PCAL 位 点 相位 误差 


Fig. 2 Phase error of single PCAL point 


以 上 海天 文 天 余 山 13m 射电 望远镜 为 例 ， 每 线程 总 带宽 为 512MH ,积分 时 间 为 10s，5MHz 间隔 下 的 信 噪 比 


为 20dB ,内 有 8 个 子 通道 ,中 心 频率 分 别 为 560MHz、592MHz、624MHz、752MHz、848MHz、912MHz、976MHz 


和 1008MHz， 表 2 给 出 各 PCAL 频率 间隔 下 的 时 延 精度 。 随 着 频率 间隔 的 减 小 ,时 延 精度 也 随 之 降低 。 因 此 ， 在 


实际 测量 中 ,通过 小 间隔 PCAL 测 出 绝对 时 延 后 ， 还 需 用 该 值 带 入 大 间隔 PCAL 的 测量 值 中 解 模 糊 度 ， 以 测 得 更 精 


确 的 链 路 时 延 。 


表 2 不 同 频率 间隔 的 时 延 精度 


Table 2 Group delay accuracy of different frequency intervals 


Interval Group delay accuracy 
500KHz 135.1ps 

1MHz 95.5ps 

2MHz 67.5ps 

5MHz 42.8ps 
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4 ”PCAL 提取 软件 
在 PCAL 信号 处 理 中 ， 需 要 对 信号 做 DFT 变换 ,以 提取 PCAL 各 位 点 的 相位 。 由 于 噪声 的 存在 ， 需 要 对 信号 


进行 积分 ， 并 通过 线性 拟 合 测量 出 各 PCAL 位 点 相位 相对 于 频率 的 斜率 ， 此 时 链 路 时 延 的 测量 值 有 : 


|. dp 
Tm = 2nf, dn 


其 中 : N 为 DFT AR , ARER , PILI DFT 点 的 斜率 

如 图 3 所 示 ， 本 文 设计 的 软件 支持 多 通道 Mark6 格式 ， 并 兼容 原始 采样 数据 RAW 格式 。 软 件 为 全 图 形 界面 ， 
操作 简单 ， 自 动 寻找 PCAL 位 点 ， 并 对 位 点 相位 进行 提取 和 线性 拟 合 。 对 于 Mark6 格式 ， 软 件 会 自动 识别 数据 由 
头 ， 根 据 帧 头 信息 自动 匹配 采样 率 和 采样 深度 ， 并 自动 分 割 通道， 用 户 只 需 按 提示 选取 文件 和 通道 。 对 于 原始 数据 
RAW 格式 ， 用 户 也 仅 需 手工 输入 采样 率 、 采 样 比特 深度 等 少量 信息 即 可 。 

要 对 绝对 时 延 进行 定 标 ， 由 于 模糊 度 的 存在 , 需要 逐 级 缩小 PCAL 频率 间隔 , 直到 链 路 时 延 测量 值 不 再 发 生变 
化 ， 即 完全 消除 模糊 度 。 在 本 软件 的 辅助 下 ， 详 细 步 骤 如 图 4 所 示 。 最 后 ， 将 5MHz 频率 间隔 下 的 链 路 时 延 测量 什 
加 上 M 倍 ( 自然 数 ) 模糊 度 ， 直 到 与 定 标的 真实 链 路 时 延 基本 相同 ， 即 可 记录 下 M 的 值 ， 该 值 为 解 算出 的 模糊 度 ， 


在 一 定时 间 内 都 不 会 发 生变 化 ， 可 以 在 后 续 观测 中 直接 使 用 。 
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3 PCAL 提取 软件 操作 界面 


Fig. 3 Ul of PCAL extraction software 
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4 绝对 时 延 定 标 步 又 


True group delay 
Solving ambiguity 


Fig. 4 Steps of finding absolute group delay 
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5 实测 
为 了 对 本 方法 进行 验证 ， 本 文选 择 上 海天 文 台 余 山 13m 射电 望远镜 的 超 宽带 数字 变频 器 ( UDC ) , 分 两 组 进行 
实测 。 其 中 一 组 实测 目标 为 对 UDC 的 绝对 时 延 进行 定 标 ; 另 一 组 依次 插入 不 同 长 度 的 线 缆 ， 模 拟 真 实 天 线 链 路 ， 
分 别 对 线 缆 和 整体 链 路 的 绝对 时 延 进行 定 标 ， 并 用 矢量 网 络 分 析 仪 对 所 定 标的 线 缆 绝 对 时 延 结 果 进 行 验证 。 
记录 器 选择 9 个 20MHz 带宽 的 子 通道 ， 每 个 子 通道 的 中 心 频率 为 : 
f.) = 500 + 50j (MHZ); j = 1,2...9 


提取 每 个 子 通道 的 PCAL 相位 后 对 所 有 子 通道 进行 带宽 综合 .如 图 5 所 示 ARAPERA (D E 2789 
整数 倍 ， 让 所 有 通道 的 PCAL 相位 基本 在 一 条 直线 上 ， 再 通过 拟 合 得 到 群 时 延 ， 并 求 出 每 个 相位 的 拟 合 残 差 。 测 试 
均 选 用 5MHz 的 PCAL 频率 间隔 所 测 数据 ,同时 按照 图 4 的 步骤 逐 级 降低 频率 间隔 来 验证 并 消除 模糊 度 。 

实测 结果 如 表 3 所 示 ， 其 中 UDC 和 全 和 链 路 绝对 时 延 仅 用 本 文 所 述 的 PCAL 定 标 方法 进行 测量 ， 依 次 插入 


的 线 线 时 延 与 采用 矢量 网 络 分析 仪 验证 的 结果 基本 一 致 ， 误差 约 为 1ns ,考虑 到 测试 环境 、 对 接头 等 差异 ， 视 为 正 
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5 带宽 综合 后 的 PCAL 相位 及 群 时 延 


Fig. 5 PCAL phase and group delay after bandwidth synthesis 


表 3 绝对 链 路 时 延 实测 结果 


Table 3 Results of absolute group delay 


Absolute delay (ns) 


PCAL VNA 
UDC 11.12 - 
全 链 路 175.02 - 
线 缆 1 14.12 15.29 
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线 缆 2 49.24 50.45 


线 绕 3 149.78 150.58 


6 结语 

本 文 提 出 一 种 用 PCAL 信号 定 标 绝对 链 路 时 延 的 方法 ,该 方法 不 需要 参考 变频 器 , 同时 用 改变 频率 间隔 的 方法 
解决 了 模糊 度 问 题 ， 可 直接 对 整个 天 线 链 路 绝对 时 延 进行 定 标 ， 用 于 满足 精准 测量 UT1、 台 站 钟 差 补偿 、 维 护 诊断 
和 提升 未 来 深 空 探测 精度 的 需求 。 经 实测 与 比 对 验证 ， 该 方法 能 准确 定 标 绝对 时 延 ， 在 10s 积分 下 达 40ps 测量 精 
度 ， 可 以 满足 需求 。 

同时 , 本文 设 计 了 一 款 用 于 台 站 快速 诊断 和 绝对 时 延 定 标的 轻 量 级 PCAL 提取 软件 ， 经 验证 ， 该 软件 能 够 解码 


Mark6 格式 和 RAW 格式 ， 并 快速 准确 提取 PCAL 相位 ， 目 前 已 在 台 站 实际 调试 工作 中 应 用 。 
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Method and application of measuring absolute link delay by PCAL 
Chang Jie'.?, Wang Jindqing134, Jiang Yongchen!, Shu Fengchun' 
(1.Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai 200030, China; 2.University of 
Chinese Academy of Sciences, Beijing 100049, China; 3. Key Laboratory of Radio Astronomy, Chinese Academy 
of Sciences, Nanjing 210008, China; 4. Shanghai Key Laoratory of Space Navigation and Positioning Techniques, 

Shanghai 200030, China) 

Abstract very long baseline interferometry (VLBI) uses phase calibration (PCAL) signal to calibrate link delay. 
However, the existing PCAL signal can only obtain the relative change of link delay. To meet the requirements of 
accurate measurement of absolute UT1, station clock error compensation, maintenance and diagnosis, and 
improving the accuracy of deep space TT&C in the future, absolute link delay is necessary to know. In actual link, 
it is difficult to measure the absolute delay due to the existence of converter and other devices. This paper 
proposes a method to measure the absolute link delay using PCAL based on the linear variation of comb 
spectrum phase with frequency. At the same time, this paper designs a lightweight PCAL extraction software for 
fast diagnosis and absolute time delay measuring, which has been applied in actual work. 


Key Words link delay, very long baseline interferometry, PCAL, phase calibration, ambiguity 
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